208 research outputs found

    Numerical solution of a non-linear conservation law applicable to the interior dynamics of partially molten planets

    Full text link
    The energy balance of a partially molten rocky planet can be expressed as a non-linear diffusion equation using mixing length theory to quantify heat transport by both convection and mixing of the melt and solid phases. In this formulation the effective or eddy diffusivity depends on the entropy gradient, ∂S/∂r\partial S/\partial r, as well as entropy. First we present a simplified model with semi-analytical solutions, highlighting the large dynamic range of ∂S/∂r\partial S/\partial r, around 12 orders of magnitude, for physically-relevant parameters. It also elucidates the thermal structure of a magma ocean during the earliest stage of crystal formation. This motivates the development of a simple, stable numerical scheme able to capture the large dynamic range of ∂S/∂r\partial S/\partial r and provide a flexible and robust method for time-integrating the energy equation. We then consider a full model including energy fluxes associated with convection, mixing, gravitational separation, and conduction that all depend on the thermophysical properties of the melt and solid phases. This model is discretised and evolved by applying the finite volume method (FVM), allowing for extended precision calculations and using ∂S/∂r\partial S/\partial r as the solution variable. The FVM is well-suited to this problem since it is naturally energy conserving, flexible, and intuitive to incorporate arbitrary non-linear fluxes that rely on lookup data. Special attention is given to the numerically challenging scenario in which crystals first form in the centre of a magma ocean. Our computational framework is immediately applicable to modelling high melt fraction phenomena in Earth and planetary science research. Furthermore, it provides a template for solving similar non-linear diffusion equations arising in other disciplines, particularly for non-linear functional forms of the diffusion coefficient

    Linking the evolution of terrestrial interiors and an early outgassed atmosphere to astrophysical observations

    Get PDF
    A terrestrial planet is molten during formation and may remain so if subject to intense insolation or tidal forces. Observations continue to favour the detection and characterisation of hot planets, potentially with large outgassed atmospheres. We aim to determine the radius of hot Earth-like planets with large outgassed atmospheres and explore differences between molten and solid silicate planets and their influence on the mass-radius relationship and transmission and emission spectra. An interior-atmosphere model, combined with static structure calculations, tracks the evolving radius of a rocky mantle that is outgassing CO2_2 and H2_2O. Synthetic emission and transmission spectra are generated for CO2_2 and H2_2O dominated atmospheres. Atmospheres dominated by CO2_2 suppress the outgassing of H2_2O to a greater extent than previously realised, as previous studies have applied an erroneous relationship between volatile mass and partial pressure. We therefore predict more H2_2O can be retained by the interior during the later stages of magma ocean crystallisation. Furthermore, formation of a lid at the surface can tie outgassing of H2_2O to the efficiency of heat transport through the lid, rather than the atmosphere's radiative timescale. Contraction of the mantle as it solidifies gives ∼5%\sim5\% radius decrease, which can partly be offset by addition of a relatively light species to the atmosphere. We conclude that a molten silicate mantle can increase the radius of a terrestrial planet by around 5%5\% compared to its solid counterpart, or equivalently account for a 13%13\% decrease in bulk density. An outgassing atmosphere can perturb the total radius according to its speciation. Atmospheres of terrestrial planets around M-stars that are dominated by CO2_2 or H2_2O can be distinguished by observing facilities with extended wavelength coverage (e.g., JWST).Comment: 19 pages, published in A&A, abstract shortene

    Timescales of chemical equilibrium between the convecting solid mantle and over- and underlying magma oceans

    Get PDF
    After accretion and formation, terrestrial planets go through at least one magma ocean episode. As the magma ocean crystallises, it creates the first layer of solid rocky mantle. Two different scenarios of magma ocean crystallisation involve that the solid mantle either (1) first appears at the core–mantle boundary and grows upwards or (2) appears at mid-mantle depth and grows in both directions. Regardless of the magma ocean freezing scenario, the composition of the solid mantle and liquid reservoirs continuously change due to fractional crystallisation. This chemical fractionation has important implications for the long-term thermo-chemical evolution of the mantle as well as its present-day dynamics and composition. In this work, we use numerical models to study convection in a solid mantle bounded at one or both boundaries by magma ocean(s) and, in particular, the related consequences for large-scale chemical fractionation. We use a parameterisation of fractional crystallisation of the magma ocean(s) and (re)melting of solid material at the interface between these reservoirs. When these crystallisation and remelting processes are taken into account, convection in the solid mantle occurs readily and is dominated by large wavelengths. Related material transfer across the mantle–magma ocean boundaries promotes chemical equilibrium and prevents extreme enrichment of the last-stage magma ocean (as would otherwise occur due to pure fractional crystallisation). The timescale of equilibration depends on the convective vigour of mantle convection and on the efficiency of material transfer between the solid mantle and magma ocean(s). For Earth, this timescale is comparable to that of magma ocean crystallisation suggested in previous studies (Lebrun et al., 2013), which may explain why the Earth's mantle is rather homogeneous in composition, as supported by geophysical constraints

    Variable stiffness McKibben muscles with hydraulic and pneumatic operating modes

    Get PDF
    McKibben muscles have been shown to have improved stiffness characteristics when operating hydraulically. However, when operating pneumatically they are compliant and so have potential for safer physical Human Robot Interaction (pHRI). This paper presents a method for rapidly switching between pneumatic and hydraulic modes of operation without the need to remove all hydraulic fluid from the actuator. A compliant and potentially safe pneumatic mode is demonstrated and compared with a much stiffer hydraulic mode. The paper also explores a combined pneumatic/hydraulic mode of operation which allows both the position of the joint and the speed at which it reacts to a disturbance force to be controlled
    • …
    corecore